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LONG-TERM  GOALS 

The  goal  of  the  project  is  to  gain  an  understanding  of  the  density  and  velocity  fields  created  by  a  surface 
forcing,  e.g.,  buoyancy  flux  and  wind  stress,  in  partially  enclosed  seas  with  hydraulically  constrained  ex¬ 
its.  Applications  to  the  Red  and  Mediterranean  Seas  are  foremost  in  our  thinking  but  the  flows  in  fjords 
and  other  embayments  of  this  type  are  also  of  interest. 

OBJECTIVES 

We  wish  to  determine  the  scaling  laws  that  relate  the  density  and  velocity  distributions  to  the  independ¬ 
ent  variables  in  the  problem.  The  latter  would  include  the  magnitude  of  the  buoyancy  flux  and  wind 
stress  and  the  geometry  of  the  sea  and  the  exit  strait.  Hopefully  these  should  then  lead  to  a  more  com¬ 
plete  understanding  of  the  prototype  and  to  suggestions  for  suitable  field  experiments. 

APPROACH 

In  order  to  understand  the  problem  outlined  above  we  have  built  a  small-scale  laboratory  model  (2.5  m 
long)  that  simulates  most  of  the  conditions  encountered  in  the  natural  system.  This  model  is  then  oper¬ 
ated  under  a  wide  range  of  operating  conditions  of  buoyancy  flux  and  model  geometry  so  as  to  gain  a 
complete  understanding  of  both  the  scaling  laws  and  the  physical  mechanisms  that  are  important  in  this 
system. 

WORK  COMPLETED 

The  effects  of  buoyancy  flux  and  a  simple  strait  geometry  have  been  completed  (see  below).  The  latter 
consisted  of  a  simple  sill  with  or  without  a  lateral  contraction.  We  have  also  run  a  number  of  tests  with  a 
contraction  alone.  We  have  also  looked  at  the  effects  of  extra,  mechanical,  mixing  in  the  main  channel  in 
order  to  try  to  understand  the  conditions  required  to  reach  the  “overmixed”  state  of  Stommel  and 
Farmer  (1953).  Also  we  started  to  look  into  strait  geometries  where  two  displaced  constraints  are  pres¬ 
ent  while  each  of  them  can  consist  of  a  sill  and/or  lateral  contraction. 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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RESULTS 


Single  Contraction  Experiments 

By  means  of  a  long  channel,  experiments  are  conducted  to  investigate  the  influence  of  the  geometry  of 
the  strait  and  the  channel  as  well  as  the  magnitude  of  the  buoyancy  flux.  Two  different  scaling  laws,  one 
by  Phillips  (1966),  and  one  by  Maxworthy  (1994,  1997a,  1997b)  are  compared  with  the  experimental 
results.  The  experiments  suggest  a  scaling  law  for  the  distribution  of  the  buoyancy  flux,  which  is  linear 
along  the  axis  of  the  channel  and  first,  was  suggested  by  Maxworthy,  rather  than  proportional  to  a 

power  of  2/3,  which  was  suggested  by  Phil¬ 
lips.  This  result  holds  for  the  experimental 
results  and  appears  to  be  valid  for  a  number 
of  natural  systems  (i.e.  Red  Sea,  Mediterra¬ 
nean  Sea).  The  experimental  facility  consists 
of  a  long  channel  that  is  some  3  m  long,  20.5 
cm  wide  and  20  cm  deep. 

The  suggested  scaling  should  read  like 
g'  =  kBo2/3xh'4/3  where  k  is  a  proportionality 
constant  depending  on  geometric  and  dy¬ 
namic  properties  of  the  flow.  Bo  is  the  buoy¬ 
ancy  forcing,  x  is  the  longitudinal  direction 
and  h  the  total  depth  of  the  water  above  the 
sill  crest.  The  factor  k  has  the  form 
k  =  F( a)W2/3Xl/3,  where  F( a)  is  a  u-shaped 
function  depending  on  the  non-dimensional  interface  depth  a  above  the  sill  crest  with  a  minimum  of 
2.52  for  a  =  0.5  and  which  goes  to  infinity  for 
a  =  0  and  1,  IT  is  the  ratio  between  the  width 
of  the  channel  Wb  and  the  width  of  the  sill  Ws, 
and  X  is  the  ratio  between  the  length  of  the 
channel  L  and  the  water  depth  h. 

Applying  this  scaling  law  to  all  of  the  experi¬ 
ments,  which  were  conducted  for  different 
values  of  W  and  varying  values  of  X,  for  each 


W,  a  very  good  agreement  is  observed. 


The  hydraulic  control  condition  is  crucial  to  this  inves¬ 
tigation.  Two  kinds  of  controls  are  considered,  one  we 
call  basic  control  and  the  other  overmixed.  For  the  ba¬ 
sic  control  case,  only  the  side  of  the  sill  towards  the 


open  end  of  the  channel  is  hydraulically  separated  from  the  flow  conditions  through  the  contraction  and 
a  <  0.5.  In  the  overmixed  case  both  sides  of  the  sill  are  hydraulically  separated  from  the  fluid  reservoirs 
and  a  =  0.5. 

In  the  figure  G"  is  the  composite  Froude  number  of  the  two  fluid  layers  as  defined  in  Lawrence  (1990) 
or  Armi  (1986),  for  example.  A  flow  region  with  Froude  number  G~  =  1  represents  a  super  critical  flow 
and  serves  to  hydraulically  separate  the  subcritical  flow  regions  (G  <  1). 

From  our  experiments  it  was  found,  that  for  basic  control  conditions  (be)  the  proportionality  constant  is 
k  =  1.04  while  for  the  overmixed  flow  (om)  k  =  2.62  W2/3X1/3. 

A  given  system  can  be  evaluated  to  whether  it  is  overmixed  or  not  by  defining  the  ratio 
r  =  g’omJg’bc  =  2.52  W2/3X1/3,  which  only  depends  on  geometric  quantities.  A  value  of  r  >  1  identifies  a 
flow  which  is  overmixed. 

To  prove  this  a  set  of  experiments  with  mechanical  mixing  at  various  locations  along  the  channel  were 
conducted. 

Single  Contraction  with  Mechanical  Mixing 

The  physical  setup  for  these  experiments  is  the  same  as  in  the  former  experiments  except  that  a  rotating 
mixer  is  used  to  force  turbulent  mixing  between  the  inflowing  and  outflowing  fluid  layers.  The  rate  of 
mixing  is  controlled  by  the  RPM  of  the  mixer  and  the  location  is  varied  x  =  50  cm  and  x  =  L. 

The  results  suggest  that  g'(L)  indeed  is  independent  of  the  specifics  of  the  mixing  within  the  channel. 
Different  rates  from  weak  to  very  strong  mixing  do  not  effect  g'(L)  which  is  especially  true  for  mixing 
locations  x  <  L  -  50  cm.  For  example  for  a  genuinely  non-overmixed  channel  the  interface  shows  strong 
turbulent  mixing  between  the  location  of  the  mixing  device  and  x  =  0  cm.  This  observation  is  backed  by 
g'(x)  ~  x  which  is  the  appropriate  scaling  for  this  case  and  this  part  of  the  channel.  Within  the  remain¬ 
der  of  the  channel  from  the  mixing  location  to  x  =  L,  g'  adjusts  so  that  g’(L)  takes  on  a  value  according 
to  the  basic  control  case. 

Moving  the  location  of  mechanical  mixing  very  close  to  the  contraction  location  (x  =  L)  generates  dis¬ 
turbances  that  effect  the  flow  conditions  across  the  contraction.  Thus  it  is  not  clear  yet  whether  it  is 
possible  at  all  to  force  a  system  that  by  default  is  basic  hydraulic  control  to  become  overmixed. 

Dual  Contraction  Experiments 

In  natural  systems  in  general  the  shallowest  and  narrowest  section  of  a  strait  do  not  coincide  as  assumed 
in  the  experiments  above.  Therefore  a  set  of  experiments  was  conducted  where  two  contraction  sections 
some  distance  apart  were  employed.  The  geometrical  properties  for  both  contractions  have  to  be  such 
that  neither  of  them  generates  overmixed  conditions  in  order  to  observe  a  feedback  between  the  two. 

The  experiments  conducted  so  far  include  a  sill  with  no  lateral  contraction  as  the  outer  most  contraction 
towards  the  open  end  of  the  channel  (i  =  1).  The  second  contraction  (1=2)  located  between  the  first 
contraction  and  the  buoyancy  sources.  This  consisted  of  different  width  ratios  W2  for  the  lateral  con¬ 
traction  while  the  sill  height  S2  was  varied  from  ,v2  =  9  to  ,v2  >  .v/. 


The  scaling  function  g'(L)  seems  to  be  according  to  the  one  found  earlier  for  the  basic  control  case  ex- 
cept  with  an  extra  factor  k  =  W2  .  The  value  of  h  to  be  chosen  is  the  smaller  one  derived  from 
h  =  H  -  si,2- 

The  underlying  physical  process  and  the  theoretical  derivation  need  yet  to  be  completely  identified. 
Theoretical  descriptions  of  width  changes  in  a  single  layer  open  channel  flow  incorporate  the  factor  of 
W  (Woodward,  1951).  The  effort  is  under  way  to  identify  the  processes  for  a  two-layer  exchange  flow 
with  width  changes  as  well. 

IMPACT/APPLICATIONS 

As  noted  above  the  laboratory  model  appears  to  support  the  claim  that  the  surface-to-depth  den¬ 
sity  difference  has  a  linear  dependence  with  distance  from  the  closed  end  of  the  sea.  This  is  consistent 
with  existing  measurements  during  some  portion  of  the  seasonal  cycle.  Also,  the  model  shows  the  im¬ 
portance  of  intense  mixing  at  the  closed  end  of  the  sea  and  the  details  of  the  hydraulic  flow  at  the  nar¬ 
row/shallow  strait.  These  observations  need  to  be  tested  during  a  field  program  that  can  use  the  experi¬ 
mental  results  to  concentrate  on  a  few,  well-chosen  locations  for  detailed  measurement  rather  than  try 
to  survey  the  whole  sea  at  low  spatial  resolution. 
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